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ABSTRACT 

We present an analysis of the scaling relations between X-ray properties and Sunyaev- 
Zel'dovich (SZ) parameters for a sample of 24 X-ray luminous galaxy clusters observed with 
Chandra and with measured SZ effect. These objects are in the redshift range 0.14-0.82 and 
have X-ray bolometric luminosity L > 10 45 erg s _1 , with at least 4000 net counts collected 
for each source. We perform a spatially resolved spectral analysis and recover the density, 
temperature T and pressure profiles of the intra-cluster medium (ICM), just relying on the 
spherical symmetry of the cluster and the hydrostatic equilibrium hypothesis. The combined 
analysis of the SZ and X-ray scaling relations is a powerful tool to investigate the physical 
properties of the clusters and their evolution in redshift, by tracing out their thermodynam- 
ical history. We observe that the correlations among X-ray quantities only are in agreement 
with previous results obtained for samples of high-z X-ray luminous galaxy clusters. On the 
relations involving SZ quantities, we obtain that they correlate with the gas temperature with 
a logarithmic slope significantly larger than the predicted value from the self-similar model. 
The measured scatter indicates, however, that the central Compton parameter yo is a proxy 
of the gas temperature at the same level of other X-ray quantities like luminosity. Our results 
on the X-ray and SZ scaling relations show a tension between the quantities more related to 
the global energy of the system (e.g. gas temperature, gravitating mass) and the indicators of 
the structure of the ICM (e.g. gas density profile, central Compton parameter yo). Indeed, by 
using a robust fitting technique, the most significant deviations from the values of the slope 
predicted from the self-similar model are measured in the L — T,L — M to t, M gas ~T, yo — T 
relations. When the slope is fixed to the self-similar value, these relations consistently show a 
negative evolution suggesting a scenario in which the ICM at higher redshift has lower both 
X-ray luminosity and pressure in the central regions than the expectations from self-similar 
model. These effects are more evident in relaxed clusters in the redshift range 0.14-0.45, 
where a more defined core is present and the assumed hypotheses on the state of the ICM are 
more reliable. 

Key words: galaxies: clusters: general - cosmic microwave background - cosmology: obser- 
vations - X-ray: galaxies: clusters 



1 INTRODUCTION 

Clusters of galaxies represent the largest virialized structures in the 
present universe, formed at relatively late times. The hierarchical 
scenario provides a picture in which the primordial density fluc- 
tuations generate proto-structures which are then subject to grav- 
itational collapse and mass accretion, producing larger and larger 
systems. The cosmic baryons fall into the gravitational potential of 
the cluster dark matter (DM) halo formed in this way, while the 
collapse and the subsequent shocks heat the intra-cluster medium 
(ICM) up to the virial temperature (0.5 < T < 10 keV). 

In the simplest scenario which neglects all non-radiative pro- 
cesses, the gravity, which has not preferred scales, is the only re- 
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sponsible for the physical properties of galaxy clusters: for this rea- 
son they are expected to maintain similar properties when rescaled 
with respect to their mass and formation epoch. This allows to 
build a very simple model to relate t he physical parameters of clus- 
ters: t he so-called self-similar model jKaiscr 1986 ; Evr ard & Henrvl 
1991). Based on that, we can derive scaling relations (see Sect.[3]l 
between X-ray quantities (like temperature T, mass A-/, entropy 
and luminosity L), and between X-ray and Sunyaev-Zel'dovich 
(SZ) measurements (like the Compton-y parameter), thanks to 
the assumption of spherical collapse for the DM halo and hy- 
drostatic equilibrium of the gas within the DM gravitational po- 
tential. These relations provide a powerful test for the adiabatic 
scenario. In particular, in the recent years the studies about the 
X-ray scaling laws (see, e.g., Allen & Fabianl I 1998; Markevitch 
ll998l : lEttori et alj2004bl ; lArnaud et al.l2005l : lvikhlinin et al.l2005t 
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iKotov & Vikhlini i32005h . together with observations of the entropy 
distribution in galaxy clusters (see, e.g JPonman et al J 19991120030 
and the analysis of simulated systems includi ng cooling and extra 
non-gravitational energy injection (see, e.g., iBorgani et al]|2004) 
have suggested that the simple adiabatic scenario is not giving an 
appropriate description of galaxy clusters. In particular the most 
significant deviations with respect to the self-similar predictions 
are: (i) a lower (by ~ 30 — 40 per cent) normalization of the 
M — T relation in r eal clusters with respect to adiabatic simulations 
fevrard etal1ll996h : (ii) steeper slopes for the M - T and L - T 
relations; (hi) an entropy ramp in th e central regions of clusters 
(see, e.g., IPonman et alJI 199*91 120031) . These deviations are likely 
the evidence of non-radiative processes, like non-gravitational heat- 
ing due to energy injection fr om supernovae, AGN, star formation 
or ga l actic winds ( see, e . g., IPearce et alj|200ll; Ifozzi & Normanl 
200ll: iBialek et alj l200ll; babul et all 120021 : IBorgani et alJ 1200a 
Brigh enti & Mathews! 120061) or cooling (see, e.g.. iBrvanl 2000). 
More recently some authors pointed out that there is a mild depen- 
dence of the X-ray scaling relations on the redshift, suggesting that 
there s hould be an evolut ion of these non-gravitational processes 
with z ( lEttori et alj|2004bl) . 

An additional and independent method to evaluate the role of 
radiative processes is the study of the scaling re lations based on 
the thermal SZ effect dSunvaev & Zeldovich|[l97()h . which offers a 
powerful tool for investigating the same physical properties of the 
ICM, being the electron component of cosmic baryons responsi- 
ble of both the X-ray emission and the SZ effect. The advantage 
of the latter on the former is the possibility of exploring clusters 
at higher redshift, because of the absence of the cosmological dim- 
ming. Moreover, since the SZ intensity depends linearly on the den- 
sity, unlike the X-ray flux, which depends on the squared density, 
with the SZ effect it is possible to obtain estimates of the physi- 
cal quantities of the sources reducing the systematic errors origi- 
nated by the presence of sub-clumps and gas in multi-phase state 
and to study in a complementary way to the X-ray analysis the ef- 
fects of extra-physics on the collapse of bary ons in cluster dark 
matter halos, both via numerical simulations ( White et al] |2002| ; 
Ida Silva et ai1l2004l ; iDiaferio et ai]|2005l ; lNagaJl2006l) and obser- 
vationally (Cooray 1999; McCarthy et al. 2003a,b; Benson et al. 
2004; LaRoque et al. 2006; Bonamente et al. 2006). 

The main purpose of this paper is to understand how these SZ 
and X-ray scaling relations evolve with redshift. In particular we 
want to quantify how much they differ from the self-similar ex- 
pectations in order to evaluate the amplitude of the effects of the 
non-gravitational processes on the physical properties of ICM. An- 
other issue we want to debate is which relations can be considered a 
robust tool to link different cluster physical quantities: this has im- 
portant consequences on the possibility of using clusters as probes 
for precision cosmology. To do that, we have assembled a sample 
of 24 galaxies clusters, for which measurements of the Compton-y 
parameter are present in the literature. Respect the previous works 
we have done our own spatially resolved X-ray analysis recovering 
X-ray and SZ quantity necessary to investigate scaling relations. 
We have performed a combined spatial and spectral analysis of the 
X-ray data, which allows us to derive the radial profile for temper- 
ature, pressure, and density in a robust way. These results, which 
have high spatial resolution, rely only on the hydrostatic equilib- 
rium hypothesis and spherical geometry of the sources. Moreover 
we can compare the observed physical quantities with the results of 
hydrodynamical numerical simulations in a consistent way. 

The paper is organized as follows. In Sect.[2]we introduce our 
cluster sample and we describe the method applied to determine 



the X-ray properties (including the data reduction procedure) and 
the corresponding SZ quantities. In Sect. [3] we report our results 
about the scaling relations here considered, including the presenta- 
tion of the adopted fitting procedure. Sect.|4]is devoted to a general 
discussion of our results, while in Sect.[5]we summarize our main 
conclusions. We leave to the appendices the discussion of some 
tecnical details of our data reduction procedure. 

Hereafter we have assumed a flat hCDM cosmology, with 
matter density parameter Qo m = 0.3, cosmological constant 
density parameter Qa = 0.7, and Hubble constant Hq = 
70km/s/Mpc. Unless otherwise stated, we estimated the errors 
at the 68.3 per cent confidence level. 



2 THE DATASET 
2.1 Data reduction 

We have considered a sample of galaxy clusters for which we have 
SZ data from the literature and X-ray data from archives (see Tables 
[T]and[2] res pectively). In particular , we have considered the original 
sample of iMcCarthv et alj J2003bh . to which we added two more 
objects from the sample discussed bv lBenson et al.1 j2004l) . For all 
these clusters we have analyzed the X-ray data extracted from the 
Chandra archive. In total we have 24 galaxy clusters with redshift 
ranging between 0.14 and 0.82, emission-weighted temperature in 
the range 6-12 keV and X-ray luminosity between 10 45 and 10 4B 
erg s _ . In the whole sample we have 1 1 cooling core clusters and 
13 no-cooling core ones (hereafter CC and NCC clusters, respec- 
tively) defined according to the criterion that their cooling time in 
the inner regions is lower than the Hubble time at the cluster red- 
shift. 

We summarize here the most relevant aspects of the X-ray data 
reduction procedure. Most of the observations have been carried 
out using ACIS-I, while for 4 clusters (A1835, A370, MS0451.6- 
0305, MSI 137.5+6625) we have data from the Back Illuminated 
S3 chip of ACIS-S. We have reprocessed the event 1 file retrieved 
from the Chandra archive with the CIAO software (version 3.2.2) 
distributed by the Chandra X-ray Observatory Centre. We have 
run the tool aciss_proces_ events to apply corrections for 
charge transfer inefficiency (for the data at 153 K), re-computation 
of the events grade and flag background events associated with col- 
lisions on the detector of cosmic rays. We have considered the gain 
file provided within CALDB (version 3.0) in this tool for the data 
in FAINT and VFAINT modes. Then we have filtered the data to 
include the standard events grades 0, 2, 3, 4 and 6 only, and there- 
fore we have filtered for the Good Time Intervals (GTIs) supplied, 
which are contained in the fltl.fits file. We checked for un- 
usual background rates through the script analyze.ltcrv, 
so we removed those points falling outside ±3<r from the mean 
value. Finally, we have applied a filter to the energy (300-9500 keV) 
and CCDs, so as to obtain an events 2 file. 



2.2 Spatial and spectral analysis 

The images have been extracted from the events 2 files in the en- 
ergy range (0.5-5.0 keV), corrected by using the exposure map to 
remove the vignetting effects, by masking out the point sources. So 
as to determine the centroid (x c ,yc) of the surface brightness we 
have fitted the images with a circular one-dimens ional (ID) isother- 
mal /3-model (fCavaliere & Fusco-Femi ano 1976), by adding a con- 
stant brightness model, and leaving x c and y c free as parameters in 
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Figure 1. The radial profiles for the projected temperature T spcc , normalized using the cooling-core corrected temperature T ow , and for density are shown for 
all objects of our sample in the left and right panels, respectively. Solid and dashed lines refer to clusters with or without a central cooling flow, respectively 



the best fit. We constructed a set of n (n ~ 15 — 40) circular an- 
nuli around the centroid of the surface brightness up to a maximum 
distance i? sp at (also reported in Table[2j, selecting the radii accord- 
ing to the following criteria: the number of net counts of photons 
from the source in the (0.5-5.0 keV) band is at least 200-1000 per 
annulus and the signal-to-noise ratio is always larger than 2. The 
background counts have been estimated from regions of the same 
exposure which are free from source emissions. 

The spectral analysis has been performed by extracting the 
source spectra from n* (n* ~ 3 — 8) circular annuli of radius r m 
around the centroid of the surface brightness. We have selected the 
radius of each annulus out to a maximum distance -R spcc (reported 
in Table O, according to the following criteria: the number of net 
counts of photons from the source in the band used for the spectral 
analysis is at least 2000 per annulus and corresponds to a fraction 
of the total counts always larger than 30 per cent. 

The background spectra have been extracted from regions of 
the same exposure in the case of the ACIS-I data, for which we 
always have some areas free from source emission. Conversely, for 
the ACIS-S data we have considered the ACIS-S3 chip only and we 
have equally used the local background, but we have checked for 
systematic errors due to possible source contamination of the back- 
ground regions. This is done considering also the ACIS "blank-sky" 
background files, which we have re-processed if their gain file does 
not match the one of the events 2 file; then we have applied the 
aspect solution files of the observation to the background dataset 
by using repro ject_events, so as to estimate the background 
for our data. We have verified that the spectra produced by the two 
methods are in good agreement, and at last we decided to show only 
the results obtained using the local background. 

All the point sources has been masked out by visual inspec- 
tion. Then we have calculated the redistribution matrix files (RMF) 
and the ancillary response files (ARF) for each annulus: in partic- 
ular we have used the tools mkacisrmf and mkrmf (for the data 



at 120 K and at 1 10 K, respectively) to calculate the RMF, and the 
tool mkarf to derive the ARF of the regions. 

For each of the n* an nuli the spectr a have been analyzed by 
using the package XSPEC ( Arnaud 1996) after grouping the pho- 
tons into bins of 20 counts per energy channel (using the task 
grppha from the FTOOLS software package) and applying the 
X 2 -statistics. The spectr a are fitted with a single-temperat ure ab- 
sorbed MEKAL model dKaastrd[l99l ; iLiedahl et al.lll995h multi- 
plied by a positive absorption edge as described in lVikhlinin et all 
(2005): this procedure takes into account a correction to the effec- 
tive area consisting in a 10 per cent decrement above 2.07 keV. The 
fit is performed in the energy range 0.6-7 keV (0.6-5 keV for the 
outermost annulus only) by fixing the redshift to the value obtained 
from optical spectroscopy and the absorbing equivalent hydrogen 
column density Nu to the value o f the Galactic neutral hyd rogen 
absorption derived from radio data dPickev & Lockmani ri990). ex- 
cept for A520, A697, A2163, MS 11 37.5+6625, MS1358.4+6245 
and A2390, where we have decided to leave TVh free due to the 
inconsistency between the tabulated radio data and the spectral fit 
result. Apart for these objects where also the Galactic absorption 
is left free, we consider three free parameters in the spectral anal- 
ysis for m— th annulus: the normalization of the thermal spectrum 
Km esc f rig dV, the emission- weighted temperature T* TO > m ; the 
metal licity Z m retrieved b y emp loying the solar abundance ratios 
from I Anders & Grevessd d 19891) . The best-fit spectral parameters 
are listed in Table [2] 

The total (cooling-core corrected) temperature T cw has been 
extracted in a circular region of radius R, with 100 kpc < R < 
R spec , centred on the symmetrical centre of the brightness distri- 
bution. In the left panel of Fig. [TJ we present for all clusters of our 
sample the projected temperature profile (T spcc ) normalized by T cw 
as a function of the distance from the centre R, given in units of 
-R2500, where R2500 is the radius corresponding to an ovedensity 
of 2500. 
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Table 2. The X-ray properties of the galaxy clusters in our sample. For each object different columns report the name, the redshift z, the identification number 
of the Chandra observation, the used ACIS mode, the exposure time t cxp , the neutral hydrogen absorption TVh (the labels / and t refer to objects for which 
Nn has been fixed to the Galactic value or thawed, respectively), the physical scale corresponding to 1 arcmin, the maximum radii used for the spatial and for 
the spectral analysis (iJ S pat and i? spC c, respectively), the emission-weighted temperature T cw , the metallicity Z (in solar units), a flag for the presence or not 
of a cooling core (labeled CC and NCC, respectively), the mass-weighted temperatu re T mw , the gas ma ss M gas , and the bol ometric X-ray lumino sity L. The 
last three columns refer to an overdensity of 2500. Sources extracted from the Mc Carthy et alj (2003b) sample and from the Benson et al. 12004) sample are 
indicated by apices (1) and (2), respectively. 
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2.3 Spectral deprojection analysis 

To measure the pressure and gravitating mass profiles in our clus- 
ters, we deproject the projected physical properties obtained with 
the spectral analysis by us ing an updated an d extended version of 
the technique presented in lEttori et al.1 J2002h and discussed in full 
detail in Appendix A. Here we summarize briefly the main charac- 
teristics of the adopted technique: (i) the electron density n e (r) is 
recovered both by deprojecting the surface brightness profile and 
the spatially resolved spectral analysis obtaining a few tens of ra- 
dial measurements; (ii) once a functional form of the DM density 
profile p — p(r, q), where q = (51, g2, ... qh) are free parameters 
of the DM analytical model, and the gas pressure Po at R spe c are 
assumed, the deprojected gas temperature, T(q, -Po), is obtained by 
integration of the hydrostatic equilibrium equation: 



P(r,q,Po) =Po 



n SBS (s)pm H 



G M(q, 



da, (1) 



where /1 = 0.6 is the average molecular weight, ran is the pro- 
ton mass. So T(q, Po) = P(q, Po)/«gas expressed in keV units. 
In the present study, to parametrize the cluster mass distribution, 
we consider two models: the universal density profile proposed by 



Navarro et al. dl997t) (hereafter NFW) and the one suggested by 
Rasia et alj j2004h (hereafter RTM). 



p(x) = 



where 



The NFW profile is given by 
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of the universe at redshift z, H z = E z Ho, E z 
+ zf + (1 - «M - Ov)(l + zf + Q A ] 1/2 , and 



Oc,NFW = 
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where c = r v i T /r s is the concentration parameter, r s is the scale 
radius, x = r/r v i r , x s = r s /r v i r . 

The RTM mass profile is given by: 

Pc,z 5c, RTM 



p(x) = 



(4) 



x(x + x*) s / 2 ' 

with i s * = rs/r v i r , where r$ is a reference radius and <5 c ,rtm is 
given by: 

5c, rtm = z — — — -r~nz — : — -ttt^ z — , ,m ■ (5) 



6 [(l + 2 3 ; 8 *)/(l+x s *)i/ 2 -2^i/ 2 ] 
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Table 1. The SZ parameters for the galaxy clusters in our sample. For 
each object different columns report the name, the central value (j/o) of the 
Compton ^-parameter, the SZ flux integrated up to an overdensity of 2500 
and over a fixed solid angle Q = 1 arcmin (j/2500 an d Vsi> respectively) 
divided by the function g(x) (see eq.[5}, and the parameter r\ (see text). For 
two objects (namely A1 914 and RX J2228+2037) the corresponding errors 
are not provided by Mc Carthy et alj i2003bh : in the following analysis we 
will assume for them a formal \u error of 20 per cent. 
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So we have q = (c, r s ) and q = (x s , r 2 oo) for the NFW and RTM 
models, respectively. 

The comparison of the observed projected temperature pro- 



Figure 2. Example of temperature spectral deprojection for cluster A1413. 
We display the two quantities which enter in the eq. lA7l in the spectral de- 
projection analysis to retrieve the physical parameters: the observed spec- 
tral projected temperature T* Q - m (stars with errorbars) and the theoretical 
projected temperature (triangules, indicated as T pro j , m in Appendix A). We 
also show the theoretical deprojected temperature T(q, Pq) (points), which 
generates T pro j im through convenient projection tecniques. 



quantities represent the 68.3 per cent confidence level and are com- 
puted by looking to the regions in the parameter space where the 
reduction of y 2 with respect to its minimum value Xmm ls smaller 
than a given threshold, fixed according to the number of degrees of 
freedom d.o.f. (see, e.g. JPress et al .1 19921) . Notice that we included 
in the eq.lQJ the statistical errors related to measurement errors of 



file Tp ro ; m (Sect. I2.2t with the deprojected T(q, Po) (eq. IA7l in 2.4 Determination of the X-ray properties 



Appendix A), once the latter has been re-projected by correcting 
for the temperature gr adient along the line of sight as suggested in 
iMazzotta et aljj2~004t) . provides the best estimate of the free param- 
eters (q, Po) through a x 2 minimization, and therefore of T(q, Po) 
(see an example in Figure[2]l. 

In the right panel of Fig. Q] we present the density pro- 
files (plotted versus r/p25oo) as determined through the previous 
method. In general, we find there is no significant effect on the 
determination of the physical parameters when adopting the two 
different DM models. Hereafter we will use the physical parame- 
ters determined using the RTM model, reported with their corre- 
sponding errors in Table [2] where we also list the exposure time, 
the number and the instrument (ACIS-I or ACIS-S) used for each 
of the Chandra observations. 

Finally we computed the total mass enclosed in a sphere of 
radius Pa as M(q)(< Pa) = f^ A p{ r , <l) dV where the ra- 
dius Pa corresponds to a given overdensity A: we considered the 
cases where the overdensity is equal to 2500 and 500. The values 
for masses and radii, together with the parameters (q, Po) for the 
RTM model, are reported in Table [3] The errors on the different 



The bolometric X-ray luminosity L(< Pa) has been calculated by 
correcting the observed luminosity 1/(100 kpc < r < P spC c) de- 
termined from the spectral analysis performed by XSPEC exclud- 
ing the central cooling region of 100 kpc (the results are reported 
in Table[2]l: 



L{< P A )=i(100kpc < r < R B 



JZ*(l+x 2 )- 3l3 x 2 dx 



(6) 



where x = r/r c , xa = Pa/t c , r c and (3 are the best-fit parame- 
ters of the /3-model on the image brightness, K COII is the normaliza- 
tion of the thermal spectrum drawn with XSPEC, and corrected for 
the emission from the spherical source up to 10 Mpc intercepted 
by the line of sight: P corr = + x 2 )~ 3/3 x 2 dx + + 

x 2 )~^x 2 (l - cos6)dx ~ J^(l + x 2 )- :i0 x 2 (l - cos0*)dx, 
with = arcsin(a;i/a;), 9* = arcsinfxo/x), xq = 100kpc/r c , 
xi = Pspcc/^c andx 2 = 10Mpc/r c . 

The gas mass M gas enclosed in a circular region having over- 
density A has been computed from the total gas density n gas j, that 
we directly obtained from the spectral deprojection, up to P spC c- 
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Table 3. Different physical properties for the clusters in our sample. For each object the different columns report the name, the minimum value for \ 2 (with 
the corresponding number of degrees of freedom d.o.f.), the virial radius r v ; r , the reference scale x B , the value of the pressure Po, the mass and the radius 
corresponding to an overdensity of 2500 (M2500 ar >d ^?2500> respectively), the mass and the radius corresponding to an overdensity of 500 (M500 and R500. 
respectively). All quantities are derived by assuming the RTM model. 
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We have checked that the exclusion of the central cooling region 
does not significantly affect the resulting values for Af gas . 

Finally we have estimated the total mass-weighted tempera- 
ture: 



(7) 



^T 3 (q,P )m ! 



which can be compared to the total emission-weighted tempera- 
ture T cw ; p represents the number of annuli inside ifesoo- Notice 
that our average deprojected temperature profile implies the follow- 
ing relation between the maximum, the deprojected and the mass- 
weighted temperatures: T max : T ow : T mw = 1 : 0.67 : 0.69 
(1 : 0.83 : 0.88 for the CC-only subsample). The physical parame- 
ters obtained in this way are also listed in Table[2]for all clusters of 
our sample. 

2.5 Determination of the Sunyaev-Zel'dovich properties 

The thermal SZ JSunvaev & Zeldovichl [l97r3) effect is a very 
small distortion of the spectrum of the cosmic microwave back- 
ground (CMB), due to the inverse Compton scatter between cold 
CMB photons and hot ICM electrons (for recent reviews see , e.g., 
lBirkinshawll999l : ICarlstrom et al.l200ztlRephaeli et al.l2005T) . This 
comptonization process statistically rises the photon energy, pro- 
ducing a distortion of the CMB black-body spectrum. The final re- 
sult is a decrease (increase) of the CMB flux at frequencies smaller 
(larger) than about 218 GHz. The amplitude of this effect is directly 
proportional to the Compton parameter y(6), which is defined as 



!/(*) = 



(Tt 



Pe(r) dl . 



(8) 



where 6 is the angular distance from the cluster centre, <jt is the 
Thomson cross-section, and P e (r) = n e (r)kT e (r) is the pressure 
of the ICM electrons at the position r; the integral is done along the 
line of sight. 

The SZ effect can be expressed as a change in the brightness: 



A/sz = 9( x ,T e ) lo y , (9) 

respectively; here Io = 2(A;Tcmb) 3 /(/ic) 2 , x = hv/kTcMB, 
Tqmb is the present CMB temperature and the function <7( x ,T e ) 
is given by: 



9(x,T e ) = 



(e x - l) 2 V e 



+ 1 



4 )(l + o(x,T e )) 



(10) 



and accounts for the frequency dependence of the SZ effect; 
the term o(x, T e ) represents the relativistic correction (see, e.g., 
Iltoh et al]|l998f) . which, however, is negligible for clusters having 
T < 10 keV. 

We consider the Compton-y parameter integrated over the en- 
tire solid angle (and given in flux units) j/a defined as: 



(11) 



va = i y(e)dn- 



To remove the dependence of j/a on the angular diameter distance 
d a (z) we use the intrinsic integrated Compton parameter Y, de- 
fined as: 



Y = d a («)j/A. 



(12) 
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The same quantity, but integrated over a fixed solid angle Q, can be 
similarly written as: 



(13) 



We fixed Q, = 1 arcmin, that is < than the field of view of OVRO, 
used i n the observations of most of the sources in our sample (see, 
e.g jMcCarfhv et alj2003ah . Notice that in order to remove the fre- 
quency dependence we have normalized Y, j/a and yn to g/ x> T e )' 
To integrate eqs. dill ) and d!3t we have recovered y(9) from 
eq. [8] by using the pressure profile P(q, Po) determined in the 
spectral analysis (Sect. [23}, renormalized in such a way that y(0) 
equals the central Comption parameter yo taken from the literature. 
This method can lead to systematics on yn and Y due to the fact 
that, even if we are assuming the true pressure profiles P(r) in 
eq. l[8j, yo has been obtained by assuming an isothermal /3-model 
inferred from the brightness profile. The value of yo is thus poten- 
tially dependent on the un derlying model of P(r). As discussed 
in rec ent works (see, e.g., [LaRoque et al. 2006; Bonamen te et al.l 
2006), the relaxation of the isothermal assumption should apply 
to the analysis of both X-ray and SZ data, to obtain a robust and 
consistent description of the physics acting inside galaxy clusters. 
Unfortunately, we have only the central Compton parameter, and 
not the complete uv— data, which are not public available: so it is 
very difficult to quantify the amplitude of this systematics, being 
yo determined through a best fit in the uv— plane. 
Nevertheless, we can give an estimate in this way: we have com- 
puted the central Compton parameter j/q iX inferred from the X-ray 
data by parametrizing first P(r) in eq. ([§) with a /3-model inferred 
on the brightness images: 



Vo 



(It 

m e c 2 



n kT gi 



(14) 



with no = n gas (r = 0) derived from the brightness profile B(r 

1 . if,,. 2\ 1/2-3/3 



B(r = 0) = 



47r(1 + , ) 4 ^A0.82ng J dx(l + x 2 ) 



(15) 



where A is the X-ray cooling function of the ICM in the cluster rest 
frame in cgs units (erg cm 3 s" 1 ) integrated over the energy range of 
the brightness images (0.5 — 5 keV). Then we have calculated j/q x 
by accounting in eq.[8]for the true pressure profile P(q, Po) recov- 
ered by the spectral deprojection analysis (Sect. l231 >. and therefore 
we determined the ratio n — yox/Vox- We notice that the pa- 
rameter n differs from the unity of < 25 per cent, comparable to 
statistical errors. 

The different quantities related to the SZ effect are listed in 
TableQ]for all clusters in our sample. 



3 THE X-RAY AND SZ SCALING RELATIONS: THEORY 
AND FITTING PROCEDURE 

3.1 The scaling relations in the self-similar model 

The self-similar model (see, e.g., Ka iseiill986h gives a simple pic- 
ture of the process of cluster formation in which the ICM physics 
is driven by the infall of cosmic baryons into the gravitational po- 
tential of the cluster DM halo. The collapse and subsequent shocks 
heat the ICM up to the virial temperature. Thanks to this model, 
which assumes that gravity is the only responsible for the observed 
values of the different physical properties of galaxy clusters, we 
have a simple way to establish theoretical analytic relations be- 
tween them. 



Numerical simulations confirm that the DM component in 
clusters of galaxies, which represents the dominant fraction of the 
mass, has a remarkably self-similar behaviour; however the bary- 
onic component does not show the same level of self-similarity. 
This picture is confirmed by X-ray observations, see for instance 
the deviation of the L — T relation in clusters, which is steeper 
than the theoretical value predicted by the previous scenario. These 
deviations from self-similarity have been interpreted as the ef- 
fects of non-gravitational heating due to radiative cooling as well 
as the energy injection fro m supernovae, AGN, star formation 



or ga l actic winds (see, e. g.. iTozzi & Normar] 



200ll: feialek et 



al. 
)4; 



| 200ll; iBorgani et al.ll2002l; Babul et alj|2002t iBorgani et al.11200- 
iBrighenti & Math ews 2006) which make the gas less centrally con- 
centrated and with a shallower profile in the external regions with 
respect the DM component. Consequently, the comparison of the 
self-similar scaling relations to observations allows us to evaluate 
the importance of the effects of the non-gravitational processes on 
the ICM physics. 

For Y and yn we have the following dependences on the cos- 
mology: 

E Z A\ /2 Y oc (p^A^yo) (e z AI /2 Ra z Y , 



and 



c -1a-1/2 77. -1a -1/2 

E z A z yn oc E z A z ' 



yo 



(16) 
(17) 

respectively, where the factor A z — 200 x 
[1 + 82 (O a - 1) / (18tt 2 ) - 39 (O z - l) 2 / (18tt 2 )] , with 
tl z — flo m (l + z ) 3 /E 2 , accounts for evo lution of clusters in an 
adiabatic scenario dBrvan & Normanlll998h . 

Assuming the spherical collapse model for the DM halo and 
the equation of hydrostatic equilibrium to describe the distribution 
of baryons into the DM potential well, in the self-similar model the 
cluster mass and temperature are related by: 

E z A] /2 M tot oc T 3/2 ; 



so we have Ea, oc (M/(p CjZ A z )) 



1/3 



(18) 

T^E^A* 1 ' 2 . By 



setting f z = E Z (A Z / A) 1 ^ 2 , from the previous e quations we 
can easily obtain the following relations (see, e.g., Markevitch 
| l998l:rA llen & F abiarJl998l:lEttori et al.l2004bl : lArnaud et alj2005l; 
piaferio et all 120051 ; IVikhlinin et al.l £0051 : iKotov & Vikhlinirj 
120051) : 



U (Y) cx (f-'yof 3 
yn oc y , 

r-l rp3/2 

fz Vo oc T 1 , 

fz 1 yo OC fgMtat , 

f7 1 yo cx (A _1 ^) 3/4 , 

fz Y oc T 5/2 , 

/ z Yoc(/ z M tot ) 5/3 , 

fzYocif-'Lf*. 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 
(26) 



We also remember here that for galaxy clusters sim- 
ilar scaling laws exist also in the X-ray band (see , e.g. , 



Ettori et alj [2004 j; I Arnaud et alj [20051 ; IKotov & Vikhlinirj l2005l 
Vikhli nin et alj2006l) : 



f-'LtxT 2 



(27) 
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f z Mt t oc T g 3 a /2 , 
f~ x L « (/ z M tot ) 4/3 . 

f A f rv T 3 / 2 
Jz Jw gas J- g as , 

f~ l L oc (/ z M gas ) 4/3 



(28) 
(29) 
(30) 
(31) 



In our work we have considered all the physical quantity at fixed 
overdensity (A z = A), i.e. f z = E z in the above equations. 



3.2 Fitting the scaling relations 

We describe here the method adopted to obtain the best-fitting pa- 
rameters in the scaling relations. Since they are power-law rela- 
tions, we carry out a log-log fit: 



log(Y) = a + Alog{X) 



(32) 



where X and Y represent the independent and dependent variables, 
respectively (hereafter Y \ X); a and A are the two free parame- 
ters to be estimated. However, in the considered scaling relations 
it is unclear which variable should be considered as (in)dependent. 
Moreover both X- and Y-data have errors due to measurement un- 
certainties, plus an intrinsic scatter. For these reasons, the ordinary 
least squares (OLS) minimization approach is not appropriate: in 
fact it does not take into account intrinsic scatter in the data, and it is 
biased when errors affect the independent variable. So we decided 
to use the BCES (Bivariate Correlated Errors and intrinsic Scatter) 
(Y\X) m odification or the bisector m odification BCES (Y, X) pro- 
posed by Akrita s~& Bershadvl dl996l) . for which the best-fit results 
correspond to the bisection of those obtained from minimizations 
in the vertical and horizontal directions. Both these methods are ro- 
bust estimators that take into account both any intrinsic scatter and 
the presence of errors on both variables. 

The results for the best-fit normalization a and slope A for 
the listed scaling relations are presented in Table [4] where we also 
report the values of the total scatter 

1/2 

S= ^(logY.-a-AlogX,- 



(33) 



and of the intrinsic scatter S calculated as: 



S = 



f(l°g Yj ~a~ AlogX,) 2 - ef ogK J / 



1/2 



(34) 



where ei og y. ■ = ey / ( Y, In 10) , with ey being the statistical error 
of the measurement Yj, and v is the number of degrees of freedom 
{y — N — 2, with TV equal to total number of data). 

Notice that in these fits the physical quantities (L, M to t , M gas , 
Y) refer to -R2500 estimated through the mass estimates based on 
the RTM model. 



3.3 On the evolution of the scaling relations 

We can extend the previous analysis by investigating the redshift 
evolution of the scaling relations at z > 0.1. Note that only two 
objects are available at 2 > 0.6 and that all CC clusters are at 
redshift below 0.45. We parametrize the evolution using a (1 + z) B 
dependence and put constraints on the value of B by considering a 
least-square minimization of the relation 



This is obtained by defining a grid of values of B and looking for 
the minimum of a x 2 -like function, defined as: 



[log Yj - a - A\ogXj - 731og(l + ^)] 2 
+ 4 log 2 X, 



log Yj 



A 2 e 2 



(36) 



log(Y) 



■A\og{X) + B\og{l + z) 



(35) 



the sum is over all data, and ei og x e ex/ (X In 10) and ei og y = 
ey/(Ylnl0) are related to the uncertainties on X and Y, re- 
spectively. The best-fit parameters values calculated by using this 
method are reported in Table [5] Again in these fits, which refer to 
same scaling relations presented in Table[4] the physical quantities 
(L, Mtot, M gaa , Y) refer to -R2500, and masses are computed by 
assuming the RTM model. 



4 DISCUSSION OF THE RESULTS 

We present here a general discussion of our results concerning the 
scaling relations. In particular we have chosen to consider both the 
whole sample (CC plus NCC objects) and the CC-only subsam- 
ple: this is done to allow a more direct comparison of our results 
with most of the works present in the literature, which are based on 
CC-sources only. Moreover this allows also to obtain at the same 
time more general relations which can be useful for future much 
extended (X-ray and SZ) cluster surveys, in which the distinction 
between relaxed and unrelaxed systems will be not easy. 



4.1 The X-ray-only scaling relations 

In this section we consider the scaling relations involving quanti- 
ties extracted from the X-ray data only. We start by examining the 
relation between Mtot and T, and finding in general a good agree- 
ment between our best-fitting slopes and the values expected in the 
self-similar model. Then we will consider the other X-ray relations, 
finding slopes which are steeper than expected from the self-similar 
model. In particular the L — T, L — M gas and Af gas — T relations 
display deviations larger than 2a, while for the L — Mtot relation 
we found agreement between the observed slope and the expected 
one. 



4.1.1 The Mtot - T relation 

Without any assumption for models on the gas density and (depro- 
jected) temperature profile, we have supposed that the DM density 
profile is well described by an analytical model (RTM or NFW). 
Thanks to the results of numerical simulations, we know, indeed, 
sufficiently well the DM physics which is in fact very simple, 
only depending on the gravity, unlike the physics of the baryons, 
which is also affected by further sources of non-gravitational en- 
ergy. Moreover we have removed the observational biases in the 
determination of the deprojected temperature (and consequently of 
the mass) by adopting the spectral-like temperature estimator (see 
Sect. 12. 3 b - In this way we have a bias-free estimate of the depro- 
jected temperature and, therefore, of the cluster mass. Below we 
focus our attention on T m w, because it is directly related to the 
total energy of the particles and so comparable to the results of hy- 
drodynamical simulation s, unlike T* e w, which is affected by obser- 
vational biases (see, e.g.. [pardini et al.l2004 : Mazzott a et alj2004l 
iMathiesen & EvrarcfeOOll) . 

First, we notice that the two different models for the DM 
profile give slightly different results. Nevertheless, at an overden- 
sity of A = 2500 the masses determined by using RTM are 
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Table 4. Best-fit parameters for the scaling relations computed by using the cluster quantities evaluated at -R2500 \ masses are estimated using the RTM profile. 
For each relation we give the logarithmic slope A (compared to the theoretically expected value A*), the normalization a, the intrinsic scatter S and the 
logarithmic scatter of the data S. The results are given both for a subsample including the CC clusters (11 objects), and for the whole sample (24 objects). In 
the column "method", symbols (1) and (2) indicate if the fit has been performed by adopting the BCES (Y \ X) or BCES (Y, X) methods, respectively. With 
the notation (yo,— 4, yn, Yg), L44, TV, M14, we indicate the Compton parameter, the X-ray luminosity, the temperature and the mass, in units of (10~ 4 , mJy, 
10 8 mJy Mpc 2 ), 10 44 erg s _1 , 7 keV, 10 14 M Q , respectively. 





Cooling core clusters 






All clusters 












11 objects 








24 objects 








relation (Y - X) 


A/ A* 


(X 


s 


S 


A/ A* 


a 


5 


S 


method 


fzYs — ft 3/0,-4 


1.22(±0.15)/1.67 


-1.07(±0.06) 


0.090 


0.113 


1.19(±0.20)/1.67 


-0.91 (±0.06) 


0.116 


0.137 


(1) 


yn - 2/0,-4 


n n^/_l_n i oo 

u.yj(±u.i4)/i.uu 


0.0l(±U.U4) 


0.033 


0.076 


0.92(±0.26)/1.00 


0.66(±0.08) 


0.120 


0.140 


(1) 


fz 2/0,-4 — T OWj 7 


2.21(±0.32)/1.50 


0.19(±0.05) 


0.138 


0.154 


2.06(±0.23)/1.50 


0.15(±0.03) 


0.123 


0.141 


(2) 


fz 1 2/0,-4 - fzM to t,14 


1.25(±0.30)/1.00 


-0.50(±0.22) 


0.248 


0.257 


1.22(±0.29)/1.00 


-0.41(±0.17) 


0.211 


0.222 


(2) 


fz 1 2/0,-4 - fz 1 L44 


0.75(±0.07)/0.75 


-0.69(±0.11) 


0.156 


0.170 


0.61(±0.05)/0.75 


-0.48(±0.07) 


0.124 


0.142 


(2) 


fz Y s - T CWi 7 


2.74(±0.23)/2.50 


-0.83(±0.03) 


0.103 


0.124 


2.64(±0.28)/2.50 


-0.74(±0.03) 


0.139 


0.157 


(2) 


/* Ys ~ fzM totM 


1.56(±0.29)/1.67 


-1.70(±0.21) 


0.235 


0.245 


1.48(±0.39)/1.67 


-1.42(±0.21) 


0.288 


0.297 


(2) 


fz Ys - fz 1 L44 


0.92(±0.11)/1.25 


-1.92(±0.16) 


0.183 


0.196 


0.81(±0.07)/0.75 


-1.58(±0.10) 


0.237 


0.248 


(2) 


fz yn — T CWj 7 


1.98(±0.46)/1.50 


0.80(±0.05) 


0.143 


0.158 


2.15(±0.45)/1.50 


0.79(±0.05) 


0.167 


0.182 


(2) 


fz 1 yn — /zM to t,i4 


1.12(±0.31)/1.00 


0.17(±0.22) 


0.239 


0.249 


1.07(±0.17)/1.00 


0.31(±0.10) 


0.278 


0.288 


(2) 


fz yn - fz L44 


0.68(±0.09)/0.75 


-0.02(±0.14) 


0.160 


0.174 


0.74(±0.10)/0.75 


0.00(±0.15) 


0.233 


0.244 


(2) 


fz 1 L44 — T cw ,7 


2.98(±0.53)/2.00 


1.18(±0.05) 


0.182 


0.183 


3.37(±0.39)/2.00 


1.03(±0.05) 


0.220 


0.221 


(2) 


I/44 — / z M totj i4 


1.71(±0.46)/1.33 


0.24(±0.32) 


0.205 


0.206 


2.03(±0.54)/1.33 


0.08(±0.32) 


0.269 


0.270 


(2) 


fzM to t,14 — T CW: 7 


1.74(±0.25)/1.50 


0.56(±0.03) 


0.000 


0.098 


1.69(±0.40)/1.50 


0.47(±0.04) 


0.044 


0.142 


(2) 


fzM to t,14 — T mW] 7 


1.63(±0.25)/1.50 


0.54(±0.04) 


0.000 


0.087 


1.69(±0.34)/1.50 


0.45(±0.03) 


0.000 


0.125 


(2) 


/z-Mgas,13 — T CWj 7 


1.94(±0.21)/1.50 


0.57(±0.02) 


0.083 


0.086 


2.09(±0.23)/1.50 


0.51(±0.02) 


0.107 


0.110 


(2) 


fz L44 — /zM gaSj i3 


1.55(±0.13)/1.33 


0.30(±0.09) 


0.083 


0.085 


1.64(±0.13)/1.33 


0.19(±0.09) 


0.131 


0.132 


(2) 



in perfect agreement with the ones determined by using NFW 



0.540 ± 0.037 and A%§ M = 1.630 ± 0.253, a£ 



0.546 ± 0.035 and A$g w = 1-590 ± 0.250). At A = 500 the 
situation becomes less clear, because for most of the clusters we 
needed to extrapolate from R spBC (corresponding to A ~ 1000) 
up to A = 500, being R spC c of order of (l/3)-(l/2) of the virial 
radius (roughly corresponding to -R2500 — fiiooo)- Hereafter we 
consider only the RTM model, even if most of the results present in 
the literature are usually based on the NFW one. 

Considering the whole sample, we find a normalization (a = 
0.45 ± .03), which i s ~ 1 (~ 5) per cent smal ler than the value 
found by I Allen et ai] feOOll) dArnaud et af] (^005)), who only con- 
sider relaxed clusters. Our normalization (a = 0.54 ± 0.04) is 
instead ~ 10 (~ 15) per cent larger than the value of Allen et al. 
(Arnaud et al.) if we only consider the CC-only subsample. This 
suggests a different behaviour depending on the presence or not 
of a cooling core (see also the left panel Fig. [3}: in fact we find 
that at A = 2500 the normalization of the NCC subsample at 
A/2500 = 5 x 1O 14 M0 (corresponding to our median value for 
the mass) is w 10 per cent smaller than for the CC-only subsam- 
ple; conversely at A = 500 the two subsamples give consistent 
normalizations, but the robustness of this result is affected by the 
fact that in this case we have to extrapolate the mass profile out of 
the region covered by observational data. 

Some other authors (e.g.. lArnaud et al]|2005l) prefer to mask 
out the central region (up to 0.1 xi?20o) in the determination of the 



mass profile. We have decided to check the effects of the inclusion 
of the cooling region in our analysis by comparing the values of the 
mass obtained by excluding or not the central 100 kpc in the deter- 
mination of the best fit parameters of the RTM profile: we pointed 
out that accounting for the cooling region does not involve any sys- 
tematic error on the determination of the mass, indeed we obtain 
more statistically robust results. 

Consequently the disagreement between CC and NCC clus- 
ters is probably due to a different state of relaxation, namely that 
the former are more regular and with more unifor m physical prop- 
erties than the latter dDe Grandi & Mol endi 2002); this is true even 
if we have masked out the most evident substructures. Notice that 
the observed mismatch is only marginally statistically significant 
(~ 1 — 1.5a). For a couple of clusters, namely A520 and A2163, 
we find that the exclusion of the unrelaxed central regions avoids 
observational biases due to the presence of local substructures: in 
particular the mass of the first (second) object increases by a fac- 
tor of ~ 2 (~ 1.5) when excluding the central 300 (360) kpc. For 
other clusters which are evidently unrelaxed, we did not find any 
convenient way to avoid possible biases: even after masking out 
the most visible substructures, the analysis of the density and de- 
projected temperature profiles still reveals the possible presence of 
local irregularities (a sort of local 'jumps' in the profiles), which 
are difficult to individuate in the brightness image. 

At A = 2500, the best fitting normalization obtained 
considering the whole sample is ~ 30 per cent below the 
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Table 5. Best-fit parameters for the redshift evolution of the scaling relations. Again, the quantities are evaluated at -R2500 an d masses are estimated by using 
the RTM profile. For each relation we list the redshift evolution parameter B, the logarithmic slope A (compared to the theoretically expected value A*), the 
normalization a, the minimum value of the function x 2 ar >d the number of degrees of freedom (d.o.f.). The results are given both for a subsample including the 
CC-only clusters (11 objects), and for the whole sample (24 objects). With the notation (3/0 —4> 3/n, Yg), L44, I7, M14, we indicate the Compton parameter, 
the X-ray luminosity, the temperature and the mass, in units of (10~ 4 , mJy, 10 8 mjy Mpc 2 ), 10 44 erg s — 1 , 7 keV, 10 14 Mq, respectively. 





Cooling core clusters 
1 1 objects 




All clusters 
24 objects 




relation (Y - X) 


B A/A* a 


Xli n (d.o.f.) 


B A/A* a 


xLn (d-o-f.) 


fz Yg-fz Vo-4 
ysi - 2/0,-4 


2.36± u ;gg 1.02(±0.09)/1.67 -1.24(±0.04) 
A - 56 -o.lo 0.85(±0.08)/1.00 0.80(±0.04) 


15.9(8) 
5.1(8) 


0.76±ojf 1.15(±0.08)/1.67 -1.00(±0.03) 
-'■ 24 ^0.24 0.82(±0.07)/1.00 0.83(±0.03) 


87.8(21) 
88.6(21) 



3/0,-4 - 



fz 1 3/0,-4 
-1 , 



- /zM tot ,14 

fz 1 3/0,-4 — fz 1 La 
fz Yg — T eWi7 
fz Yg - f z M totM 

fz Yg — fz 1 L 4 4 



fz 

fz 
t-1 

fz yn 



yn 
yn 



■ fzM to t,14 

■ fz 1 L44 



-2.12 
-2.44 
0.04 
-1.08 
-2.32 

2.40; 

-4.00 
-4.88 
-1.52 



+0.96 
0.96 
1.08 
2.52 
+0.48 
0.48 
-1.12 
-1.16 
-1.96 
-2.60 
+0.44 
0.48 
+ 0.96 
0.96 
1.68 
2.44 
+0.48 
0.44 



2.41(±0 
1.35(±0 
0.69(±0 
2.98(±0 
1.68(±0 
0.70(±0 
2.31(±0 
1.33(±0 
0.58(±0 



25)/1.50 
23)/1.00 
05)/0.75 
31)/2.50 

25) /1.67 
05)/1.25 

26) /1.50 
22)/1.00 
05)/0.75 



0.40(±0.03) 
-0.33(±0.19) 
-0.57(±0.07) 
-0.71 (±0.04) 
-1.55(±0.20) 
-1.81(±0.07) 
1.20(±0.03) 
0.52(±0.18) 
0.34(±0.07) 



23.4(8) 
29.6(8) 
48.0(8) 
10.6(8) 
20.5(8) 
58.9(8) 
16.6(8) 
26.9(8) 
47.5(8) 



+0.36 
0.32 
+0.48 
0.48 
+0.16 
-0.16 
0.40 
0.40 
+0.52 
0.52 
0.20 
0.16 
+0.36 
0.32 
+0.44 
0.40 
, -,,+0.16 
1 -' z -0.16 



0.08 
0.08 
-0.32 
0.28 
0.28 
0.12 
-1.44 
-1.52 



2.08(±0 
0.98(±0 
0.62(±0 
2.66(±0 
1.14(±0 
0.62(±0 
1.95(±0 
0.82(±0 
0.40(±0 



17) /1.50 
10)/1.00 

03) /0.75 
20)/2.50 
13)/1.67 

04) /0.75 

18) /1.50 
12)/1.00 
03)/0.75 



0.12(±0.02) 
-0.27(±0.07) 
-0.44(±0.05) 
-0.78(±0.03) 
-1.26(±0.08) 
-1.35(±0.05) 
0.97(±0.02) 
0.64(±0.08) 
0.65(±0.05) 



82.4(21) 
55.6(21) 
99.7(21) 
37.9(21) 
78.9(21) 
206.0(21) 
71.2(21) 
120.0(21) 
178.0(21) 



fz L44 — T CWi 7 

fz L44 — / z M tot ,i4 

fzM to t,14 — Tew, 7 

/zMtot.14 — r mw ,7 

fzM gaSi i3 — T cw .7 
fz L44 — / z M gaSi i3 



H>.96 
-1.00 
H.24 
-1.56 
4+1.12 
-1.20 
0+1.24 
— 1 32 

0.16+ ' 56 



-O.72: 
-1.32j 
0.56^" 
-0.8 



-0.92 



0.60 
+0.24 
0.24 



3.27(±0.29)/2.00 
1.29(±0.16)/1.33 
1.79(±0.30)/1.50 
2.00(±0.28)/1.50 
2.00(±0.16)/1.50 
1.43(±0.03)/1.33 



1.25(±0.03) 
0.66(±0.13) 
-1.00(±0.28) 
-1.09(±0.27) 
0.57(±0.02) 
0.45(±0.03) 



69.40(8) 
5.4(8) 
11.1(8) 
7.1(8) 
39.2(8) 
73.7(8) 



,+0.52 
-0.52 
•1+0.56 
-0.56 
a+0.52 
-0.52 
-,+0.48 
-0.48 

0.84+g;^ 

-,+0.12 
-0.12 



0.92: 

-0.24^ 
-0.08^ 
-0.32^ 



-0.60^ 



4.05(±0.24)/2.00 
1.36(±0.11)/1.33 
2.30(±0.24)/1.50 
2.32(±0.22)/1.50 
2.17(±0.12)/1.50 
1.63(±0.02)/1.33 



0.85(±0.03) 
0.53(±0.07) 
-1.51(±0.22) 
-1.54(±0.21) 
0.41(±0.02) 
0.26(±0.02) 



190.0(21) 
40.5(21) 
48.4(21) 
33.0(21) 
127.0(21) 
358.0(21) 



value found in the non-ra diative hydrodynamic simulations by 
iMathiesen & Evrar3 j2001)Fl for the CC-only subsample, the nor- 
malization is ~ 20 per cent below the theoretical value. The dis- 
crepancy is slightly reduced (~ 15 — 20 p er cent) with respec t to 
the adiabatic hydrodynamic simulations bv lEvrardetaIlJl99d) . 

The picture emerging from numerical simulations with a 
more sophisti c ated I CM modeling is different. The simulation by 
Borga niet al.l J2004I) . which includes radiative processes, super- 
nova feedback, galactic winds and star formation, suggests a nor- 
malization which is in rough agreement with our whole sample, and 
15 per cent lower with respect to the CC-only subsample. Notice, 
however, that the re-analysis of the same simulation data made by 
iRasia et alj d2005l) . who adopted a different definition of tempera- 
ture, the spectroscopic-like one (which is not consistent with our 
definition of mass-weighted temperature; see above for a more de- 
tailed discussion), gives a higher (~ 40 — 50 per cent) normaliza- 
tion. 

Finally we notice that the slope of the M — T relation is, in- 
deed, in agreement with the theoretical expectations (A* — 1.5). 

Considering the results at an overdensity of 500, we found 
a good agreement (at la level) between observed and theoretical 
slopes. 



Our analysis suggests no evolution (B — —0.88 



1 We have rescaled their results from A = 500 to A = 2500. 



+ 1.24 
1.32' 



B = —0.32 ± 0.48) , in agreement with the literature (see , e.g., 
iFinoguenov et alj200ll ; lEttori et alj2004bl : I Allen et al.ll200l[) . 

We compare also our intrinsic sc atter, which is con sistent 
with zero, with the one estimated by IRasia et al", I l l2005h : they 
find a scatter of « 30(16) per cent by considering the emission- 
weighted (spectroscopic-like) temperature. We reach similar con- 
clusions compari ng our intrinsic scatter with the value retrieved by 
iMotl et all d2005h . 



4.1.2 The L - T relation 

We find (see the upper-right panel of Fig.O a marginal agreement 
of our results on t he slope of this rela tion (^4 al1 = 3.37±0.39), with 
those obtained bv lEttori et"al] j2002l) . who found A = 2.64 ± 0.64 
at A = 2500: however, their sample contains colder objects, for 
which a flatter relation wo uld be expected. O ur results also agree 
with the analysis made by Mar kevitchl Jl998t) : A = 2.64 ± 0.27. 
Notice that his cluster sample is not directly comparable with ours, 
since it covers different ranges in redshift and temperature. 

We compare our results abo ut the scatter (S = 0.220 and 
S = 0.221) with those obtained bv lMarkevitchl (l998), who found 
a smaller value: S = 0.103 (see, however, the previous comments 
on the different characteristics of the two samples). 

Moreover, we find (at ~ la) a positive (negative) redshift 
evolution for all clusters (CC-only subsample), i.e. we notice a 
mildly different behaviour on the evolution CC and NCC clusters. 
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For comparison lEttori et al.l J2002h found B 
their sample of clusters at higher redshift. 



-1.04 ± 0.32 for 



Hints of negative evolution are observed (B 



-1.32 



+ 1.24 



B n 



-0.24 ±0.56). 



Regarding the normalization we observe a slightly different 
behaviour when the CC-only subsample and whole sample are con- 
sidered: q cc = 1.18 ± 0.05 and a al1 = 1.03 ± 0.05, respectively. 
We notice that the luminosity of the CC clusters is systematically 
larger than that of the NCC clusters, even if we have corr ected it 
for the cooling flow (see Sect. l2.4t . as already observed by Fabian 
dl994l) . On the contrary numerical simulations predict that the re- 
moval of the gas from the X-r ay emitting phase reduces th e lumi- 
nosity dMuanwong et al.ll2002h . This confirms that cooling dBrvanl 
2000; Voit & Br vanl l200lh is not effective in removing baryons 
from the X-ray phase, be cause of the presence of an extra-source of 
feedback or pre-heating dBalogh et al II 19991 : ICavaliere et alii 19981 ; 
iTozzi & Normanl 1200 ll ; iBabul et al. 2002), which maintains the 
ICM at warm temperature dBorgani et al.ll2002l) . Alternatively, the 
more evident negative evolution of the CC clusters compared to the 
NCC ones (especially in the yn — X-ray(SZ) relations) could indi- 
cate different states of relaxation, bei ng the former more regular , 
relaxed and virialized than the latter dPe Grandi & Molendill2002l) . 



4. 1.3 Other X-ray scaling relations 

Here we discuss our results for the relations not shown in the fig- 
ures. For the M gas — T relation we find a ~ la discrepancy be- 
tween the slope of this relation in the CC-only subsample (^4°° = 
1.94 ± 0.21 and A M = 2.09 ± 0.23) and the theoretical expecta- 
tion for the self-similar model (A = 1.5). Nevertheless, our esti- 
mate is consistent with the results already present i n the literature. 
By ap plying a /3-model to recover the gas mass, IVikhlinin et al.1 
( 1998) measured A = 1.71 ± 0.23 at the baryon overdensity 1000 
(approximately corresponding to the virial DM overdensity). Our 
slope is also in good agreem ent with the value (A — 1.98 ± 0.18) 
found bv lMohr et all |l999), always by applying the /3-model. We 
ha ve also a marginal a greement (at la level) with the value found 
bv lEttorietai] ( l2004bl) (A = 2.37 ±0.24), wh o make use of the /3 - 
model and apply the correction for E z . Finallv lEttori et all 120021) . 
combining a spectral analysis and the application of a /3-model to 
the brightness distribution and without correcting for E z , found 
A = 1.91 ±0.29 for A = 2500 and A = 1.74±0.22 at A = 500. 
The results of this last paper also suggest a low intrinsic scatter, 
in good agreement with our analysis (S = 0.079). We point out 
here that we find some differences between CC and NCC clusters 
at A = 2500, because of the contribution of the cooling core region 
(< 100 kpc); at A = 500 this effect becomes negligible because 
the behaviour of the gas mass is dominated by the contribution from 
the external regions (Af gas oc r). Finally no significant evolution is 
observed (B = 0.16±g for the CC clusters; when we consider 
the whole sample, we notice a more significant positive evolution 
(B = 0.84 ±0.28) 

Regarding the L — M to t the best-fit slope for the CC-only 
subsample (A cc = 1.71 ±0.46) is in good agreement with the re- 
sults obtain e d bvlReiprich & B6hringeifd2002l) (A = 1.80 ±0.08) , 
lEttori et alj d2002h (A = 1.84 ± 0.23) and lEttori et al.l d2004bh 
(A = 1.88 ± 0.42). The observed scatter we measure (S cc = 
0.206, S al1 = 0.270) is slightly sm aller than in previous anal- 
ysis by iReiprich & Bohringeil d2002l) (S = 0.32), and in agree- 
ment with lEttori et alj d2002l) (S = 0.26). This seems to suggest 
that the methods we applied to correct the observed luminosity 
(see Sect. I2.4t and to determine the total mass are quite robust. 



For the L — A/ gas law we measure a slope which is discrepant 
with respect to the theoretical value expected in the self-similar 
model. This relation, together with the one between Af gas — T and 
A/tot — T, has the lowest intrinsic scatter between the X-ray only 
scaling laws. Moreover we have a significant evidence of a negative 
redshift evolution. 



4.2 The scaling relations involving the SZ effect 

In this section, we discuss first the Y — yo and yn — yo relations, 
which are linking the SZ properties only (see Fig.0, and then the 
relations between SZ and X-ray quantities (see Fig.[5]l. The impor- 
tance of these relations relies on the possibility of providing new 
insights into the general physical properties of the ICM, in a way 
complementary to the X-ray view. In particular, the different depen- 
dence on the gas density and temperature of the SZ flux (~ n e T) 
with respect to the X-ray brightness (~ n 2 T 1 ' 2 ) can allow to re- 
duce some of the biases present in the X-ray analysis. The pres- 
ence of substructures and inhomogeneities in the ICM can indeed 
strongly affect some of the X-ray determined physical parameters, 
like temperature and luminosity. An independent approach through 
the SZ analysis of some physical quantities can shed more light on 
the limits of validity of the ICM self-similar scenario. 



4.2.1 The Y — yo, yn — yo relations 

For both relations, we find slopes which are smaller than the ex- 
pected o nes. The discrepancy w e measure is larger than the one 
found bv lMcCarthv et aljd2003bl) . This is likely due to the fact that 
the self-similar model predicts a pressure profile which is steeper 
than the observed one: including extra-gravitational energy draws a 
picture in which the gas density (and consequently the pressure) has 
a profile shallower than the DM density. This is also confirmed by 
the observation that there are differences between CC (which are 
obviously more subject to non-gravitational processes) and NCC 
clusters, having the former a slightly (~ la) smaller integrated 
Compton parameter. We point out that the dispersion in these re- 
lations is very high, probably because of the systematics on the 
reconstruction of the integrated Compton parameter (see Sect. l2.5t . 

We measure a a strong negative evolutio n in the yn — yo 
relation. As pointed as iMcCarthv et al.l d2003bl) . this different be- 
haviour of the yn — yo relation (more in general of the yn— X- 
ray and yn— SZ relations) concerning the evolution is likely due 
to the fact the SZ effect within a fixed angular size samples larger 
physical region at higher redshifts. This means that the effect of 
non-gravitational processes are relatively more pronounced if the 
SZ flux is measured within smaller physical radii, where the den- 
sity of the ICM is higher: this is expected in a scenario of either 
preheating, where we can assign a fixed extra-energy per particle, 
or cooling, where the radiative cooling is more prominent in the 
denser central regions. This is also in agreem ent with the general 
picture emerging by studying entropy profiles (Ponman et al. 2003; 
IPratt et al.l 120061; I Voit & Ponniar] 120031 ; ITozzi & Normanl 1200 ll) . 
which is affected just in the central regions by non-gravitational 
processes, while the self-similarity is roughly preserved in the halo 
outskirt, where the dynamics is still dominated by the gravity. 
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Figure 3. The relations between Mtot-T mw (left panel) and L-T cw (right panel). In each panel the filled circles represent cooling core (CC) sources, while 
the diamonds are the no-cooling core (NCC) ones. The solid line refers to the best-fit relation obtained when considering all clusters of our sample, the dashed 
one represents the best-fit when the CC sources only have been considered and the dot-dashed is the best-fit obtained by fixing the slope to the self-similar 
value. 
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Figure 4. As in Fig.[5]but for the relations between Y — yo (left panel) and yn — yo (right panel). 



4.2.2 The yo —T,Y—T,yQ—T relations 

We note that yo — T is the only scaling relation that deviates by 
> 3cr from the self-similar slope (see Table 4 and 5) both when 
only CC clusters and CC+NCC objects are considered. Moreover, 
we measure an higher normalization in the CC-only subsample, 
probably due to the inclusion of the cooling regions during the SZ 
data reduction and the subsequent fit in the visibility plane. These 
results , in good agreement with the ones presented in lBenson et al.l 
d2004h . are consistently obtained with both a robust BCES fit and 
a x 2 -minimization. By applying the former technique, this rela- 
tion is also the one that shows the smaller scatter (both total and 



intrinsic) around the best-fit. Furthermore, the \ -approach indi- 
cates a significant negative evolution among the 11 CC clusters 
{B cc _ _ 2 .i2l |g at 2 - 5cr ; Xmin = 23.4 with 8 d.o.f.) that 
disappears when the whole sample of 24 objects is considered. 
For the NCC sources we do observe hints of positive evolution 
(B ncc = 0.641°;*°): this points to a different behaviour of the 
cool core and non-cool clusters in the central regions, and dif- 
ferent state of relaxation of the gas as suggested by the compar- 
ison of the normalization of the fit (a cc = 0.19 ± 0.15 and 
q ncc _ o.i4 ±0.35). 

The best-fitting relations for Y — T show a value for the slope 
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Figure 5. As in Fig.[5]but for the relations between yo — T ew (left panel), yo — Mtot (central panel), yo — L (right panel). 



in agreement with the value predicted by the self-similar model 
either when we consider the CC-only clusters or the whole sample, 
unlike for the yo — T relation: this is probably due to the sensitivity 
of yo to the cooling region. On the contrary for the yn—T relation, 
when we consider the CC clusters, we observe a good agreement 
with the self-similar predictions (A cc = 1.98 ±0.46 A* = 1.50). 

Our results confirms that the Y — T relation exhibits a smaller 
scatter than the yo — T one, as naively expected. Finally we find that 
the yet — T relation has a la rger scatter than t he Y — T one, in con- 
trast with what obtained by McCarthy et a 1. (2003a). Moreover we 
notice in the CC-only subsample a mildly larger scatter compared 
to the whole cluster sample. 

4.2.3 The yo - M, Y — M, yn — M relations 

These relations show a very good agreement between observed and 
self-similar slopes, with a scatter a factor of 2 larger than the corre- 
lation with T (see the previous subsection). We do not confirm the 
low scatter, S ~ 10—15 per cent, for t he Y — M re latio n suggested 
from the numerical simulations by [Nagai ( 200of) and iMotl et all 
j2005l) : this indicate possible bias in the determination of Y. But it 
is possible that the present simulations are not completely adequate 
to reproduce the observed quantities, being the ICM modeling in 
hydrodynamical codes quite complex. 

The normalization of the Y — Mtot relation has been in- 
vestigated in dedicated hydrodynamical simulation s to discrimi - 
nate between different ICM physics. For example, iNagail tOOdt) 
uses non-radiative (NR) and with gas cooling and star formation 
(CSF) simulated clusters to find a normalization that varies by 
about 70 per cent: for a typical cluster with M2500 = 5 x 1O 14 M0, 
Y NR = J1.32±°;J°) x 10- 4 and Y CSF = (9.0llJ;S) x 10" 5 
at z = Qj. At the same mass and overdensity, and fixing the slope 
to the self-similar model, our observed normalization is: Y cc = 
(5.32 ± 1.06) x 10~ 5 and Y al1 = (8.06 ± 1.35) x 10~ 5 for CC- 
only and all clusters, respectively. At A = 200, the observed nor- 
malizations are Y cc = (1.30 ±0.74) x 10" 5 and Y al1 = (1.22 ± 
0.53) x 10 , systematically lower than the results in Nagai (2006) 
( Y NR = 5.13lg;g x 10- 5 and Y CSF = 3.95±g-fl x 10" 5 ) 
and more in agreement with the results by Ida Silva et alj J2004h . 

2 Here we are following his definition of Y, corresponding to To = 1 in 
eq. il Q , and we adopt his cosmological parameters. 



that measure Y NR = 1.85 x 10~ 5 , Y c ° o1 = 1.73 x 10~ 5 and 
ypro-hcat _ 2 g Q x 1Q -6 for non _ ra( iiative, cooling (cool) and 
pre-heating (pre-heat) simulations, respectively. 

We obtain, therefore, that our CC clusters, for which we obtain 
the most robust estimates of the total mass at the overdensity of 
2500 (see Subsect. fc.l.lt . well reproduce the distribution measured 
in the Y — Mtot plane of the objects simulated including extra 
physical processes. Similar conclusions can be drawn for Y — T mw 
and Y — L relations. 

Finally, we find a negative evolution for the relations under ex- 
aminations at ;> lcr confidence level for the CC-only clusters (see 
Table 5). The slopes of the correlations tend, however, to deviate 
from the self-similar predictions more significantly than the mea- 
surements obtained with the robust fitting technique. If we fix the 
slope to the self-similar value A* in these relations between SZ and 
X-ray quantities, we still obtain a neg ative evolu tion at ~ 1 — 2a 
confidence level. We note here that lNagail f2006), on the contrary, 
does not find any hint of evolution in the Y — M relation. 



4.2.4 The yo — L, Y — L, yn — L relations 

In general we find a good agreement between the best-fitted slope 
and the self-similar prediction. Compared to other scaling relations, 
in these cases the intrinsic scatter is very small (~ 0.15 for the yo — 
L relation estimated in the CC-only subsample). We do not observe 
significant differences between CC and NCC clusters, being the 
estimates of luminosity corrected for the cooling core. 

Regarding the evolution, we find suggestions (at 3<r level) for 
a negative evolution in the yn — L relation (B cc — —1.52^1^). 
We observe instead positive evolution in the Y — L relation, B cc = 
2.40jlo'48> b ut negative evolution when we consider the NCC clus- 
ters (B NCC = -0.80l J 4 ). 



5 CONCLUSIONS 

We have presented an analysis of X-ray and SZ scaling relations 
of a sample of 24 galaxy clusters in the redshift range 0.14-0.82, 
selected by having their SZ measurements available in literature. 
We have analyzed the Chandra exposures for these X-ray luminous 
objects. We have reconstructed their gas density, temperature and 
pressure profiles in a robust way. Then, we have investigated the 
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scaling relations holding between X-ray and SZ quantities. By as- 
suming an adiabatic self-similar model, we have corrected the ob- 
served quantities by the factor f z = E z , neglecting the factor A z , 
checking that the final results do not change significantly in this 
way: so we can compare our results with the work in the literature. 
We have estimated the values of normalization, slope, observed and 
intrinsic scatters, and evolution to quantify the amplitude of the ef- 
fects of the non-gravitational processes in the ICM physics. In this 
sense, the combined study of the SZ and X-ray scaling relations 
and their evolution in redshift is a powerful tool to investigate the 
thermodynamical history in galaxy clusters. Indeed, the departures 
from the self-similar predictions observed in some of the scaling 
laws studied in our work confirm that the simple adiabatic sce- 
nario is not wholly adequate to describe the physics of the X-ray 
luminous clusters, because it does not account for a further non- 
gravitational energy besides the potential one. We remind that our 
results are, by construction, more robust at -R2500, where no ex- 
trapolation is required and the determinations of the mass (at least 
for CC clusters) and the reconstruction of the integrated Compton 
parameter are reliable. These results can be here summarized as 
follows. 



• We observe a good agreement of the normalization of the 
Mtot — T relation between our results and the ones obtained in 
hydrodynamical numerical simulations. The other X-ray scaling re- 
lations involving a direct propagation of the absolute value of the 
measured gas density show a steeper slope than expected from self- 
similar predictions. Departures larger than 2a are observed in the 
L-T (^ a11 = 3.37 ± 0.39 vs. A* = 2), L - M gas (A al1 = 
1.64 ± 0.13 vs. A* = 1.33) and M gas - T (A al1 = 2.09 ± 0.23 
vs. A* = 1.5) relations. These results are consistent with previ- 
ous analysis on high-z X-ray luminous galaxy clusters (see, e.g . , 



lEttori et alJ|2002l ; lKotov & Vikhlininll2005l:lMaughan et alj|2 006). 

• Correlations between the investigated SZ quantities and the 
gas temperature have the largest deviations from the slope predicted 
from the self-similar model and the lowest scatter among similar 
relations with different X-ray quantities. The measured scatter is 
comparable to what is observed in the relations between X-ray pa- 
rameters. The Y — T relation shows the lowest total and intrinsic 
scatter both when CC clusters only and the whole sample are con- 
sidered. 

• We observe a strong negative evolution in the yn —X-ray and 
yn — SZ relations. A plausible explanation is that the SZ effect 
within a fixed angular size samples larger physical region at higher 
redshifts. That means the effect of non-gravitational processes are 
relatively more pronounced within smaller physical radii. 

• The observed normalization of the Y — Mtot relation in 
cooling-core clusters at A = 2500, that provide the most robust 
estimates of the total masses in our cluster sample, agrees well 
with the predicted value from n umerical simulations (see, e.g., 
Ida Silva et al.1 |2004 iNagail l2006h . In particular, we confirm the 
trend that lower normalization are expected when some feedback 
processes take place in the cluster cores: for a cluster with typical 
M2500 » 5x 1O 14 M , we measure y cc = (5.32 ±1.06) x 10" 5 
in the sample of CC objects where the cooling activity is expected 
to be very effective, and F al1 = (8.06 ± 1.35) x 10" 5 in the whole 
sample. However, we have to note that the normalization in hydro- 
dynamical simulations is strictly related to the adopted recipes to 
describe physical processes, like gas cooling and star formation. 
These processes are also responsible for the production of the cold 
baryon fraction, the amount of which is still under debate when 



comp ared to the observational constraints (see, e.g. jBorgani et al.l 
12006b . 

• The SZ - X-ray relations are, in general, well described by a 
self-similar model parametrized through the dependence upon f z , 
when a robust fitting technique, that considers both the intrinsic 
scatter and the errors on the two variables, is adopted. On the con- 
trary, when an evolution in the form (1 + z) B is investigated by 
a x 2 -minimization with error propagations on both X and Y vari- 
ables, we measure a strong negative evolution at > la level of con- 
fidence for all relations that involve SZ quantities (yo,Y,yn) and 
the X-ray measured gas temperature and total mass. The slopes of 
the correlation tend, however, to deviate from the self-similar pre- 
dictions more significantly than the measurements obtained with 
the robust fitting technique. If we fix the slope to the self-similar 
value A* in these relations between SZ and X-ray quantities, we 
obtain stronger hints of negative evolution for the yo — Mtot re- 
lation (B cc = -0.88 ± 0.94) and for the Y - Mtot relation 



-2.30 ± 1.13). 



Our results on the X-ray and SZ scaling relations show a ten- 
sion between the quantities more related to the global energy of 
the system (e.g. gas temperature, gravitating mass) and the indica- 
tors of the ICM structure (e.g. gas density profile, central Comp- 
ton parameter yo). Indeed, by using a robust fitting technique, the 
most significant deviations from the values of the slope predicted 
from the self-similar model are measured in the L — T, L — Mtot, 
Mg as — T, yo — T relations. When the slope is fixed to the self- 
similar value, these relations show consistently a negative evolution 
suggesting a scenario in which the ICM at higher redshift has lower 
both X-ray luminosity and pressure in the central regions than the 
self-similar expectations. These effects are more evident in relaxed 
CC clusters in the redshift range 0.14-0.45, where a more defined 
core is present and the assumed hypotheses on the state of the ICM 
are more reliable. 

A likely explanation is that we need an increase in the cen- 
tral entropy to spread the distribution of the gas on larger scales: 
this could be achieved either by episodes of non-gravitational 



heating due to supernovae a nd AGN (see, e.g ., Evrard & Henrvi 
199 ll: ICavaliere et al.l Il999l: iTozzi & Norman[|200ll: iBialek et al l 



200l |: lBrighenti & Mafhewsll200ll : lBabul et alj2002l : lBorgani et al.l 

2002), or b y selective removal of low-entropy gas through coolin g 



(see, e.g. JPearce et al.l200ll : rVoit & BrvarfcOOWWu & Xuel2002h . 
possibly regulated by some mechanism supplying energy feedback 
[e.g. the semi-analytical approach proposed by IVoit et aD [2002) 
and t h e numerical simulations discussed by Muanwon g et al.l 
d2002l) : lTomatore et al.l j2003ri ; lKav et all fe(X)3h l. 
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APPENDIX A: SPECTRAL DEPROJECTION 
TECHNIQUE 

The deprojection technique decomposes the observed X-ray emis- 
sion of the i-th annulus into the contributions from the volume frac- 
tion of the j-th spherical shells with j < i, by fixing the spectrum 
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normalization of the outermost shell to the corresponding observed 
values. We can construct an upper triangular matrix , where the 
column vectors V 1 , V 2 , ...V n represent the "effective" volumes, i.e. 
the volume of the j-th shell contained inside the j-th annulus (with 
j ^ i) and corrected by the gradient of n 2 inside the j-th shell (see 
Appendix B for more detail), so as: 



Ki 



<j dV 



v#n; 



(Al) 



In the previous equation n e = (n e ,i, n e> 2, • n e ,n). being n 
the total number of annuli, having internal (external) radius 
n n ,i , n n ,2 , ■•• ,rm,n (r ou t,i , r ou t,2 , r ut,n), with n ~ 15 - 
40; Ki is the MEKAL normalization of the spectrum in the i- 
th annulus; the operator # indicates the matrix product (rows by 
columns). Notice that the integral f >i n 2 jdV is of the order of the 
emission measure inside the i-th ring0The inversion of this matrix 
allows us to determine n e ,j. 

The values of Ki are obtained by rescaling by the observed 
number of counts in the i-th ring the faked Chandra spectrum with 
absorption, temperature and metallicity measured in that ring. The 
errors are computed by performing 100 Monte Carlo simulations 
of the observed counts. We pointed out that the uncertainties in 
the estimates of the projected temperature do not reflect into high 
systematic errors in the determination of Ki, because of the mild 
dependence on T of the cooling function A(T) integrated in the 
considered band (0.5 - 5 keV): A(T) oc T~ a , with 0.1 < a < 0.2 
for T ~ 7 - 12 keV. 

This approach is very powerful, because does not require any 
"real" spectral analysis, which could suffer of the poorness of the 
statistics and would need at least ~ 2000 net counts per annulus: we 
can determine the projected density in annuli even with very small 
counts (~ 200 — 1000). In other words we have an improvement 
(of about one order of magnitude) of the spatial resolution in the 
spectral analysis. 

Concerning the temperature analysis, we have determined its 
value Tj in the j-th shell, by assuming analytic relations for the 
mass density profiles: p = p(r, q), where q = (qx, q%, ... qn) are 
suitable parameters. As discussed in Section 2.4, we consider two 
functional forms, a NFW profile with q = (c, r s ) and a RTM pro- 
file with q = (x*, r v i r )- 

We performed a spectral deprojection of the observed 
temperature T s * hcll in a set of n* annuli with width much 
larger than the previous ones, with internal (external) radius 
r* n ,i , n*n,2 , -,r*n tn * (Tovit,! , »"out,2 , -,4t,i,«) corresponding 
to the ones of the rings in which we have estimated the projected 
temperature (see Sect. 12. 2t . with n* <C n (n* ~ 3 — 8), so as 
to have at least 2000 counts per annulus. The deprojection method 
works in this way: 



-L ring 



V*# T* sll eu-e" 



/ L ring,m , 



(A2) 



where the operator "•" indicates the product: T* sric n ■ e* = 
(T*sheii,i e*, T* s hcii,2 £2, ••■T 1 *shcii,n* e*.). In eq.l|A2), e* = 



V 2 ,...V n , L* rin 



V* 1 #L r ing is the emissivity, V* = jV 1 
is the luminosity of the m-th ringQ and the generic parameter V* 

3 Hereafter we assume that the index j (i) indicates the shell (ring) of the 
source of radius (r- ln , r ou t ) ■ 

4 Hereafter we assume that the index k (m) indicates the shell (ring) having 



has the same meaning as above, but it is evaluated in n* annuli. The 
inversion of the matrix in eq. dA2t allows us to finally estimate the 
deprojected temperature T£. 

We computed the theoretical temperature Tj by numerically 
integrating the equation of the hydrostatic equilibrium (eq. [TJ, as- 
suming spherical geometry (r = r). Then we constructed a grid 
of values for Po and for the parameters q entering the DM den- 
sity profiles, so as T, = Tj(q, Po). In particular for Po we have 
considered the range Po — 3ap ^ Po ^Po + 3<rp , where Po is 
the expectation value of Po. So we can estimate the temperature 
*Tj(q,fl>) =P(r)/n gas (r). 

Since the temperature Tj(q, Po) obtained in this way is given 
on a set of n annuli with spatial resolution much better than the 
deprojected temperature T* hcll k defined in the n* annuli only, we 
perform a (gas mass-weighted) average to calculate the temperature 
T fe ave (q, P ) in the fc-th shell: 



P fe aVe (q,Po) 



r k^ r j< r k + l 



W j T j (q,P )dV J 



WjdVi 



(A3) 



where Wj — nj and dVj represents the volume of the j-th shell, so 
as to reproduce a mass-weighted temperature. A % 2 -minimization 
between Tfe Ve (q, P ) and T* hcll k (with error a T , ), 



2 I 2 

I T 



(A4) 



provides us the best estimate of (q, Po). 

We also con sidered an alternative approach to determine 
(q, Po). Following iMazzotta et alj ( f2004h . we perform a weighted 
average of Tfc(q, Po) to compute a projected spectral-like temper- 
ature T pro j im (q, Po) in the m-th ring to be compared with the ob- 
served temperature T* j m of the m-th ring: 

T proj , m (q,Po)^V*#^ vc (q,P )-^(q,Po)^ /Alng.rn, (A5) 



where Wi 



"(q,-Po) 



3/4, £ ring (q,Po) 



V*#w(q,Po); V* = [v 1 , V 2 , ...V n * , V n * +1 , V n * +h ] is an 
extension of the volume matrix V* which takes into account the 
contributions (up to a distance of 10 Mpc) coming from the h an- 
nuli external to P spat . We have to use the followi ng fitting func- 
tion, w hich is a simplified of the functional form of lVikhlinin et al.l 
j2005h : 



»(r) = 



(r/r c ^ a (l + ryr a -')- Eh 



(A6) 



(1 +r 2/ r 2)3/2/3- a /2 

with 7 = 3, and (TJ to extrapolate n e (R), the pressure and temper- 
ature in regions outside P sp at- Notice that the previous definition of 
temperature is a very powerful way to remove observational biases: 
in fact we are weighting different regions along the line of sight 
using different temperatures which are obtained by performing a 
spectral fit of a single-temperature model. With this approach we 
have a robust determination of the deprojected temperature profile. 
The best estimate of (q, Po) is obtained through a x 2 -minimization 
between T proj , m (q, P ) and T proj m : 



C^proj,m(q; Po) ^Dro^m) 



m — 1 
.2 



(A7) 



Here or proj m accounts for the statistical errors in eq.{TJ coming 
from the measured errors for n ga , B (r). The reduced \ 2 resulting 
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from this method is better than in previous case: this is likely due 
to the fact that the deprojected temperature T s * hcll strongly relies 
on assumptions, like spherical symmetry and uniform density pro- 
file, which are not completely satisfied in real clusters. Moreover 
the values of T s * hcll k are not independent: in fact we relate the de- 
convolved temperature, gas density and spectra normalization of 
the outermost shell to its observed values and then we compute the 
physical parameters in the m-th annulus by opportunely accounting 
for the contributions of the fc-th shell (k ^ m): this could propa- 
gate possible systematic errors from the external regions, where the 
determination of the physical properties cannot be so adequate be- 
cause of the bad statistic. All the deprojected quantities presented 
in the present work refer to the second approach (T pr oj ) only. 



a good spatial resolution, for example when we consider n annuli 
(n ~ 15 — 40) in the brightness image (see Section [2~3t . But this 
is false when we have n* annuli, with n* <C n (n* ~ 3 — 8) in 
the spectral analysis, for which a larger statistics (at least ~ 2000 
net counts per annulus) is required. In this last case, for example, 
it is possible to underestimate the true density in the external re- 
gions by 5 — 10 per cent by using the geometrical volume instead 
of the effective one: this corresponds to set P equal to X (being 

X the identity matrix) in eq. dB2t . and P* oc K (see eq. lAlt . The 
analysis we performed shows that the case in which adopting the 
effective volume is important is in eq. dA2t : using an effective vol- 
ume can avoid to introduce systematic errors in the determination 
of the cluster masses. 



APPENDIX B: DETERMINING THE EFFECTIVE 
VOLUME 

iKriss et ail dl983l) computed the geometrical volume of the j-th 
shell intercepted by the i-th annulus (with j i) as: 

/Tout, r(rf a . -r 2 ) 1 ' 2 

V?=4lT drr ' dz. (Bl) 

Jr in . h-rl^.-r 2 ) 1 ' 2 

Notice that when we use a geometrical volume to deproject 
the physical parameters (like the density as example) we are as- 
suming that they are nearly constant in the shell. This introduces 
a systematic bias in the deprojected quantity, that tends to be in- 
creased when the gradient of the physical parameter is not neg- 
ligible or when the rings are wide. McLaughlin (1999) partially 
corrected this bias by referring the density to an average radius, 

^((rl 2 +rL /2 )/2) 2/3 . 

Here we introduce a new definition of the volume, the effective 
volume V, which takes into account the real gradient of the physical 
parameters as a function of the radius. We assumed that we are 
weighting the unknown physical parameter V in the j-th shell using 
a function w(R), whose gradient is only due to the case of the 
squared density (w(R) oc n 2 ). 

We have modeled the density inside the j-th shell 
as a local power-law, n(R) — n ej f(R)~ a , where 
f(R) = CR/iWj), r refi = (r inj + r outj )/2, a(R) = 
-log(n J+1 /n J )/lc-g(r rcfj+1 /r rofj ) + O(o). We first calcu- 
lated a(R) by relying on the initial density obtained from the 
geometric volume-deprojection on a radius r re f, : in this way the 
introduced errors on a are negligible (0(a)). 

We define r as the projection of R on the sky plane, with R 2 = 
r 2 + z 2 , being z the distance along the line of the sight. So, if 
n(r le £. ) is the density in the j-th shell, the observed parameter V* 
is related to the theoretical one by: 

P* = j dVPw(R) = J dVp n 2 (R) = 

dVf(R)- 2a )# (tC/p) = V#(rC/p) . (B2) 
So we can re-write the effective volume V as: 
V(=[ dVf(Ry 2a =4TT [drr [ ^ * dz f(Ry 2a . (B3) 



The effective volume Vf is equal to the geometric one VI if 
a = 0, i.e. when we have negligible gradients of n(R) 2 in the j-th 
shell. This is approximately true only in the case in which we have 



